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ABSTRACT 

This paper describes the response of the IceCube neutrino telescope located at the geographic South Pole to outbursts of MeV neutrinos from the 
core collapse of nearby massive stars. IceCube was completed in December 2010 forming a lattice of 5160 photomultiplier tubes that monitor a 
volume of ~ 1 km 3 in the deep Antarctic ice for particle induced photons. The telescope was designed to detect neutrinos with energies greater 
than 100 GeV. Owing to subfreezing ice temperatures, the photomultiplier dark noise rates are particularly low. Hence IceCube can also detect 
large numbers of MeV neutrinos by observing a collective rise in all photomultiplier rates on top of the dark noise. With 2 ms timing resolution, 
IceCube can detect subtle features in the temporal development of the supernova neutrino burst. For a supernova at the galactic center, its sensitivity 
matches that of a background-free megaton-scale supernova search experiment. The sensitivity decreases to 20 standard deviations at the galactic 
edge (30 kpc) and 6 standard deviations at the Large Magellanic Cloud (50 kpc). IceCube is sending triggers from potential supernovae to the 
Supernova Early Warning System. The sensitivity to neutrino properties such as the neutrino hierarchy is discussed, as well as the possibility to 
detect the neutronization burst, a short outbreak of v c 's released by electron capture on protons soon after collapse. Tantalizing signatures, such as 
the formation of a quark star or a black hole as well as the characteristics of shock waves, are investigated to illustrate IceCube's capability for 
supernova detection. 
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1. Introduction 

On February 23, 1987, a burst of mainly electron anti- 
neutrinos with energies of a few tens of MeV emitted by 
the Sup ernova SN1987A was record ed simultaneously b y the 
Baksan (lAlekseev et all 1 19871). 1MB dBionta et all Il987l) . and 
Kamiokande-II dHirata et all 1 19871 1 19881) detectors, a few hours 
before its optical counterpart was discovered. With just 24 neu- 
trinos collected, stringent limits on the mass of the v e , its life- 
time, its magneti c moment and the n umber of leptonic flavors 
could be derived dKotake et alll2006l) . As of now, SN1987A re- 
mains the only source of neutrinos that has been detected outside 
of our solar system. Although the optical detection of supernova 
explosions has a long history, detailed features of the gravita- 
tional collapse can only be studied with neutrinos, which carry 
away nearly 99 % of the gravitational binding energy soon after 
the collapse. The current generation of detectors is capable of 
detecting many orders of magnitude more neutrinos and thus it 
can study details of the gravitational collapse and neutrino prop- 
erties. 

The rate of galactic stellar collapses, including those ob- 
scured in the optical , is estimated by various methods to be 
~ (1 - 1)1 100 years (iDiehl etall 2006; St romlll994l) . A com- 
pilation in iGiunti & Kiml d2007l) narrows the expected range to 
(1.7 - 2. 5)/ 100 years by taking into account experimental and 
theore tical limits. The bes t exper imental upper limit is < 9.3 /l 00 
years dNovoseltseva et al.ll2009l) . 

While differences in the onset of the neutrino emissio n be- 
tween various models are small (Ka chelriess et all [2005), the 
models have yet to overcome problems with the supernova ex- 
plosion mechanisms. The theoretical knowledge about the neu- 
trino emissi on at times longer than several 100 us aft er the delep- 
tonization dBuras et all 120031; iKitaura et all 120061) is limited. 
However, three characteristic phases are expected: a rapid lumi- 
nosity increase during collapse with the appearance of a shock 
breakout burst, an accretion phase ending after (9(0.5) s during 
which the neutrino flux of all flavors is maximal, and a cool- 
ing phase. The (9(20) s duration of supernova neutrino emission 
is determined by the neutrino diffusion time scale in the dense 
matter inside the proto-neutron star. The exact features will de- 
pend on the progenitor mass with modulations introduced by the 
dynamics of the collapse. 

IceCube is primarily designed to observe TeV neutrino 
sources with a wide lattice of light sensors embedded in highly 
transparent glacier ice use d as Cherenkov medium. However, 
it was recognized early by iPrvor et all d!988l) and lHalzen et al.l 
(1996) that neutrino telescopes offer the possibility to monitor 
our Galaxy for supernovae. In spite of the much lower neutrino 
energies of (9( 1 MeV) involved in a supernova burst, Cherenkov 
light induced by neutrino interactions will increase the count rate 
of all light sensors above their average value. Although the in- 
crease in the noise rate in each light sensor is not statistically 
significant, the effect will be clearly seen once the rise is consid- 
ered collectively over many sensors. Low photomultiplier noise 
rates, low photon absorption in the Cherenkov medium and a 
large number of sensors are essential. 

IceCube is uniquely suited for this measurement due to its 
location and 1 km 3 size. The noise rates in IceCube's photomul- 
tiplier tubes average around 540 Hz since they are surrounded 
by inert and cold ice with depth dependent temperatures rang- 
ing from -43 °C to -20 °C. At depths between (1450 - 2450) m 
they are also largely shielded from cosmic rays. The noise rate is 
further reduced by the use of detector components with reduced 
radioactivity. The detected signal rate is essentially independent 



of the photon scattering length and depends linearly on the ab- 
sorption length of a; 100 m in ice. 

The expected signal significance in IceCube is somewhat re- 
duced due to two types of correlations between pulses that in- 
troduce supra-Poissonian fluctuations. The first correlation in- 
volves a single photomultiplier tube. It comes about because a 
radioactive decay in the pressure sphere can produce a burst of 
photons lasting several //s. The second correlation arises from 
the cosmic-ray muon background; a single cosmic ray shower 
can produce a bundle of muons which is seen by hundreds of 
optical modules. 

The 5 160 photomultipliers are sufficiently far apart such that 
the probability to detect light from a single interaction in more 
than one DOM is small. Effectively, each DOM independently 
monitors several cubic-meters of ice. The detection principle 
was demons trated with th e AMANDA experiment, IceCube's 
predecessor dAhrens et al.Ll2002h . 

The inverse beta process v e + p — > e + + n dominates super- 
nova neutrino interactions with (9(10 MeV) energy in ice or wa- 
ter, leading to charged particle tracks of about 0.5 cm ■ E v j MeV 
length. Considering the approximate E\ dependence of the cross 
section, the light yield per neutrino roughly scales with E\. 
Due to the low rate of galactic supernovae, it is imperative that 
the detector operates stably for a long time. IceCube was de- 
signed to operate for at least 10 years and is well suited for 
such a purpose owing to an automated online data acquisition, 
analysis software and alert system. As neutrinos may escape 
from an exploding supernova at much higher matter densities 
than photons, neutrinos will be observable several hours before 
their optical counterpart. The detailed observation of the onset 
of a supernova explosion is of much interest to astronomers. 
Since 2009, IceCube has been sending real-time datag rams to the 
Super nova Early Warning System (SNEWS) dAntonioli et all 
2004) when detecting supernova candidate events. SNEWS 
has been set up to broadcast a reliable alert to the astron- 
omy community when a supernova has been detected by sev- 
eral neutrino detectors within seconds of each other. Cu rrently , 
Super -Kamiokande dFukuda et all 120021; llkeda etall 12 007), 



LVD ( Ag lietta et aTll2002l), Bor exino dAlimonti et alll2009l) and 



IceCube (Ahrens et all l2004h contribute to SNEWS, with a 
number of other neutrino and gravitational wave detectors plan- 
ning to join in the near future. 

This paper describes the technical details and expected 
physics capability of IceCube as a detector for core collapse su- 
pernovae. It also summarizes the performance of the detector 
while it was still under construction. The outline of the paper 
is as follows: Sect. [2] describes physics processes in supernovae 
for selected models and oscillation scenarios that are relevant to 
the performance studies presented in this paper, Sect.[3]describes 
the aspects of the IceCube detector relevant to the detection of 
MeV supernova neutrinos. In Sect. |4l we discuss the processes 
that lead to a detectable signal in IceCube as well as the online 
analysis that processes and monitors the data, triggers events and 
sends out alerts to SNEWS. Sect. [5] describes the performance 
of the detector over two years and the systematic uncertainties 
expected when assessing the sensitivity of the detector. Sect. [6] 
discusses IceCube's potential in the study of astrophysical and 
neutrino properties, and finally, conclusions are given in Sect. [7] 

2. Supernovae and Neutrinos 

After the core of an aging massive star ceases generating en- 
ergy and the corresponding radiative pressure from nuclear fu- 
sion processes, it undergoes a sudden gravitational collapse as 



IceCube Collaboration: R. Abbasi et al.: IceCube Sensitivity for Low-Energy Neutrinos from Nearby Supernovae 



soon as its inactive core grows beyond the Chandrasekhar mass 
limit. After several steps to relieve thermal and degeneracy pres- 
sure from the dense electron gas, the collapse stops once nuclear 
densities are reached and an incompressible proto-neutron star is 
formed. Matter falling on its surface is promptly stopped and its 
momentum is inverted forming an outward moving shock wave. 

Neutrinos of different flavors are initially trapped in their 
relative neutrino spheres as the mean free path of neutrinos is 
smaller than the size of the supernova core at densities larger 
than 10 13 kg/m 3 . The shock wave following the collapse dis- 
sociates nuclei, which suddenly increases the number of pro- 
tons, resulting in an increase in electron capture and the produc- 
tion of a burst of v e . The timescale of this neutronization burst 
("deleptonization peak") is on the order of 10 ms during which 
much of the energy driving the shock wave is carried away. The 
shock stalls but is presumably soon revived by interactions of 
the large flux of neutrinos generated in the proto-neutron star. 
The models describing the prompt neutronization burst appear 
to be robust and consistent dKachelriess et aUl2005l) . The proto- 
neutron star subsequently cools over ~ 20 s. The neutrino flux 
decreases until neutrinos are no longer produced in the cooled 
down proto-neutron star (see e.g. Fisc her et alj (|2010) and refer- 
ences therein). 

The released gravitational potential is carried away by huge 
numbers of neutrinos and to a small extent by heating and ex- 
pelling the star's outer layers. Less than 1 % of the gravita- 
tional binding energy of a supernova is emitted as kinetic en- 
ergy of matter and optically visible radiation. The remaining 
99 % is released as neutrino energy, of which about 1 % will 
be carried by electron neutrinos from the initial neutronization 
burst. Most neutrinos and antineutrinos, distributed among all 
flavors, are created during the subsequent cooling processes. An 
estimated £ tota i = J L^jt'W s 3 x 10 53 erg = 1.87 x 

10 59 MeV dBurrows et all [l992) is carried away by the in- 
tense neutrino burst produced predo minantl y through thermal 
Kelvin-Helmholtz cooling reactions dSuzukiLfl99lt) . Here Li?' al 
is the time dependent all flavor s upernova neutrino and anti- 
neutrino luminosi ty. According to Thomp son et alj d2003l) and 
Bura s et alj d2003l) . the mean energy is expected to be about (13 
- 14) MeV for v e , (14 - 16) MeV for i> e an d (20-21) MeV for 
all other flavors (v x ); the Garching model ( Kita ura et all [2006) 
differs in that the mean energies for v e and v x turn out to be ap- 
proximately equal. For a supernova at d — 10 kpc distance and 
an average neutrino energy E Y = 15 MeV, the summed flux of all 
neutrino and antineutrino types, ^SSi = E Xot . i \l(Ajid 2 E v ), flow- 
ing through the detector is * 10 16 irr 2 . Mo re on the theory of 
core collapse supernova can e.g. be found in lJanka et al.l (2006) 
and references therein. 

In the following paragraphs, we b riefly introduce the 
Law rence-Livermore ( Totani et all 1 19971) and Garching mod- 
els dKitaura et all 12006b that are used as benchmarks. In ad- 
ditio n, we introduce specifi c models by Dasgu pta et al.l (2010) 
and ISumivoshi et al.l d2007l) that were selected to demonstrate 
IceCube's physics performance in Sect. [6] We also discuss the 
effect of neutrino oscillations on the expected signals, and intro- 
duce the parametrization of the energy spectra chosen for this 
paper. 

The spherically symmetric Lawrence-Livermore simulation 
was performed from the onset of the collapse to 18 s after the 
core bounce, encompassing the complete accretion phase and a 
large part of the cooling phase. It is modeled after SN 1987 A 
and assumes a 20 M progenitor. The total emitted energy is 



2.9 x 10 53 erg, of which 16% is carried by v e with 15.3 MeV 
energy on average. 

The newer spherically symmetric Garching simulations in- 
clude more detailed information on neutrino energy spectra amd 
use a sophisticated neutrino transport mechanism. They cover 
0.80 s following the collapse of an O-Ne-Mg 8-10 M Q pro- 
genitor star, that is destabilized due to rapid electron capture 
on neon and magnesium. This class of stars may represent up 
to 3 % of all core collapse supernovae. Recen t simulations 
by dFischer et all 1201 Ot iHiidepohl etaiT, 1201 Ol) extend over 
22 s from the collapse of a 8.8 M Q progenitor to the com- 
pleted formation of the deleptonized neutron star. They are the 
only examples so far where one-dimensional simulations obtain 
neutrino-powered supernova explosions and two-dimensional 
simulations yield only minor dynamical and energetic modifca- 
tions. In Table [4] we also refer to a tw o-dimensional axisymmet- 
ric simulation bv lMarek et al.l d2009l) of a 15 M Q progenitor star 
that covers 0.38 s following the collapse. Results with full multi- 
an gle neutrino transpo rt in two dimensions have been reported 
bv lBrandt et aild2011l) . 

Following or iginal work by Takahara & Satol d!988l) . re- 
cent simulations (Dasgupta e t all 



2010) of certain stellar core- 



collapse supernovae predict a sharp burst of v e several hundred 
milliseconds after the prompt v e neutronization burst associated 
with a quark-hadron phase transition at high baryon densities. A 
detection of this prominent feature would constitute direct evi- 
dence of quark matter. 

The gravitational collapse of less than solar metallicity stars 
exceeding 25 solar masses will lead to a limited stellar explo- 
sion, while stars exceeding 40 solar masses are not expected to 
explode at all. In both cases a black hole will develop 0(1 s) af- 
ter bounce. At this point, the neutrino emission quickly comes 
to an end, providing a uniq ue signature for black hole forma- 
tion dSumivoshi et all 120071) and a mode l independent time-o f- 
flight measurement of the neutrino mass dBeacom et al. [ l200l . 

Neutrinos streaming out of the core will encounter matter 
densities ranging from 10 13 kg/m 3 to zero. Assuming an en- 
ergy of 15 MeV, they pass through an MSW-resonance layer 
at s2 x 10 6 kg/m 3 associated with Am^ or A«j 3 , depending 
on the neutrino mass hierarchy, followed by a second layer at 
« 2 x 10 4 kg/m 3 associated to the quadratic neutrino mass differ- 
ence Amj 2 . Both mix the initial fluxes of v e , v e and v x depending 
on the survival probabilities. Although the survival probabilities 
depend on the details of the density profiles and generic predic- 
tions are impossible, we consider three limiting cases as bench- 
marks to discuss the effect of the assumed neutrino hierarchy 
on the spectra observed with IceCube. Scenario A describes the 
normal neutrino hierarchy case and Scenario B represents the 
inverted hierarchy case with a static density profile of the super- 
nova, both paired with a relatively large mixing angle 813 > 0.9°. 
In Scenario C, the mixing angle #13 is assumed to be very small 
(#13 < 0.09°) and the hierarchy may be either normal or in- 
verted. One should be aware that the predictions are affected by 
the unknown density profile of the collapsing star. In addition, 
forwards or backwards running single or multiple shock waves 
or bubbles can form within the supernova, causing steep density 
gradi ents and - in some cases - change s in th e oscillation behav- 
ior d^misitiLl HOOS IChoubixitSlllQOi. 

As the neutrinos propagate through the Earth, they un- 
derg o matter induced oscillations (MSW ef fect). The neutrino 
flux dGiunti & KimLl2007HDighe et alll2004l) decreases by up to 
8 %. The effect depends sensitively on the zenith angle, the su- 
pernova model and the assumed neutrino properties. Given the 
systematic uncertainties, it will be difficult to establish this effect 
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with IceCube. We will therefore include earth oscillation effects 
in the systematic uncertainty. 

For this paper, the supernova is considered to be close to 
blackbody source of neutrinos while it is cooling down. The 
neutrino energies then follow a modified Fermi-Dirac distribu- 
tion. Many model predictions discussed in this paper adopt the 
following parametrization for the neutrino differential flux 
at the position of the Earth at distance d from the supernova: 



dE v 



where 



^SnW frJ7 -FT , 

—f(Ey, Ey, a v ) 

4nd 2 E v 



(1) 



f(E v , Ey, a v ) = I — I X r7r - — (2) 



E v {t) 



T(l + a v {t)) 



is the norma lized energy dist ribution depending on a shape pa- 
rameter a v dKeil et all 120031) . Theory provides the time de- 
pendent supernova luminosity L v SN (t) for the neutrino species 
v = v e , v e , v x with corresponding energies E v , Other model pre- 
dictions are transferred to this framework by fitting the provided 
spectra. 

3. The IceCube Detector 



IceCube (lAhrens et all [2004) was installed in the Antarctic ice 
sheet at the geographic South Pole between January, 2005 and 
December, 2010 by lowering cable assemblies, called strings, 
into holes drilled in the ice using hot water. IceCube instru- 
ments a volume of about 1 km 3 of clear ice between depths of 
1450 m and 2450 m below the surface with a coarse lattice of 
5160 Digital Optical Modules (DOMs). Each DOM consists of a 
photomultiplier tube housed in a pressure-resistant glass sphere. 
Once the water in the holes refreezes, the DOMs are embedded 
into the ice sheet with good optical coupling. The DOMs are 
installed on 86 cable strings, of which 80 will be separated by 
roughly 125 m forming a triangular grid, with each string con- 
taining 60 DOMs vertic ally spaced by 17 m. Mo re details on the 
detector can be found in Achterber g et alj (|2006). The remaining 
6 strings constitute the denser DeepCore sub-array. DeepCore 
strings are separated by approximately 60 m and are located near 
the center of IceCube. Each of these strings contains 60 high 
quantum efficiency DOMs, with the bottom 50 DOMs vertically 
spaced by 7 m and located between depths of (2107 - 2450) m 
below a dusty ice layer with reduced transparency. The remain- 
ing 10 DOMs, vertically separated by 10 m and located above 
the dust layer, instrument a volume between (1750 - 1860) m. 

IceCube is complemented by a surface array called IceTop, 
consisting of a pair of DOMs encased in ice tanks near the top of 
each string. Due to their highe r noise and sen sitivity to the fluctu- 
ating solar particle flux ( Abbasi et all [2008), the IceTop DOMs 
do not contribute to the IceCube supernova detection system. 

The DOM is the fundamental element in the IceCube archi- 
tecture. Each DOM is housed in a 13" (33 cm) diameter, 0.5" 
(1.27 cm) thick borosilicate glass pressure sphere. It contains 
a Hamamatsu R7081-02 (R7081-MOD in case of DeepCore) 
10" ( 25.4 cm) hemispherical photomultiplier tube (Abbasi et al 



2010) as well as several electronics boards containing a proces- 
sor, memory, flash file system and realtime operating system that 
allows each DOM t o operate as a compl ete and autonomous data 
acquisition system (Ab basi et all 120091) . It stores the digitized 
data internally and transmits the information to a surface data 
acquisition system on request. 



The supernova detection relies on continuous measurements 
of photomultiplier rates. The rate information is stored and 
buffered on each DOM in a 4-bit counter in 1.6384ms time bins 
(2 16 c ycles of the 40 MHz clock). The main data acquisition sys- 
tem (Abbasi et al., 2009) transfers the data asynchronously to 
the independently operating supernova data acquisition system 
(SNDAQ). For the real-time processing, the information is syn- 
chronized by a GPS clock and regrouped in 2 ms bins. 

The South Pole is out of reach for most communication 
satellites and high bandwidth connectivity is available only 
for abou t 6 hours per da y. Therefore, a dedicated Iridium- 
satellite dPratt et all 1 19991) connection is used by the SNDAQ 
host system to transmit urgent alerts. In that case, a short data- 
gram is sent to the northern hemisphere. The receiving system 
parses the message and forwards information on the supernova 
candidate event to the international SNEWS group. The time de- 
lay between photons hitting the optical module and the arrival of 
the datagram at SNEWS stands at about 6 min, providing close 
to real-time monitoring and triggering. In order to test the signal 
path, an internal trigger threshold is adjusted to transmit 1-2 
background triggers per day. 

Due to satellite bandwidth constraints, the data are re-binned 
in 0.5 s intervals and then subjected to a statistical online anal- 
ysis described in Sect. 14.21 the fine time information in 2 ms in- 
tervals is transmitted for a period starting 30 s before and ending 
60 s after a trigger flagging a candidate supernova explosion. 
The system is surveyed by the IceCube experiment control and 
monitoring system ("IceCube Live"); supernova alerts are im- 
mediately distributed by e-mail to notify experts. 

The optical and noise properties of the DOMs are crucial 
for the understanding of IceCube's supernova detection and will 
hence be discussed in more detail in the following subsections. 



3. 1 . Optical Properties of the Digital Optical Module 

The photomultiplier was chosen on the basis of low dark counts 
and good time and charge resolution. Its bialkali photocathode 
has a spectral response in the range 300 nm to 600 nm with a 
peak quantum efficiency of (25 ± 1) % at 420 nm, well-matched 
to the Cherenkov signal spectrum and the optical properties of 
the glacial ice. Dark count rates for standard efficiency DOMs 
of around 540 Hz are typical for DOMs at ice temperatures be- 
tween -43 °C at -1450 m depth and -20 °C at -2450 m depth. 
The quantum efficiency of high efficiency DOMs is roughly 1.35 
times higher, while the noise increases approximately by a fac- 
tor 1.25. The glass of the pressure sphere was selected for high 
transmission in the sensitive region of the photomultiplier and 
a low rate of background photons from intrinsic radioactivity in 
the glass. Optical transparency extends well into the near-UV, 
with 50% transmission at Tso% ~ 340 nm, but drops to a few 
percent at 310nm. 

The photomultiplier is mechanically attached to the glass 
pressure housing with a silicone elastomer gel (GE6156 RTV). 
This gel matches the refractive index of the glass (n=1.48) to 
reduce optical losses at the medium interfaces. The spectral 
transparency of the gel extends to approximately 250 nm with 
T$o% ~ 300 nm. The combined response of the glass, gel, and 
photomultiplier is a critical input to the IceCube Monte Carlo 
simulation package. The detection efficiency of the DOM to a 
head-on parallel light beam is shown in Fig.Q] 

Most of the PMTs are operated at a gain of 10 7 , so single 
photoelectrons produce pulses of approximately 8 mV amplitude 
and 10 ns width across the load impedance. The programmable 
front end pulse discriminator is set to 2mV, a factor of * 10 
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Fig. 1. Overall DOM efficiency versus wavelength for head-on 
illumination of the 0.0856 m 2 DOM cross section. The average 
value in the 300 - 600 nm range, weighted by the wavelength 
dependence of Cherenkov light emission, is « 7.1%. 
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Fig. 2. Probability density distribution of time differences be- 
tween pulses for noise (bold line) and the exponential expec- 
tation for a Poissonian process fitted in the range 15 ms < Ar < 
50 ms (thin line). The excess is due to bursts of correlated hits, 
as can be seen from the 50 ms long snapshot of hit times shown 
in the insert. 



above the RMS noise level dAbbasi et al.L l2010h . On average, 
85 % of all single photo-electron pulses pass the discriminator 
threshold. 



3.2. Noise Properties of the Digital Optical Module 

Several effects contribute to the prevailing average rate of 540 Hz 
for standard efficiency DOMs. Atmospheric muons (16 Hz), 
thermal emission from the photocathode (<10Hz), and photo- 
multiplier induced afterpulses (» 30 Hz) all play a role, but the 
majority of hits are due to radioactive decays of which a large 
fraction initiates bursts of hits lasting for up to 15 ms. These 
bursts are presumably scintillation of residual cerium in the glass 
of the photomultiplier and the pressure sphere energized by (3 
and a decays. The 40 K content in the IceCube pressure spheres 
is specified to produce less than 100 Bq leaving trace elements 
from uranium and thorium decay chains as the main source of 
radioactivity. 

The main characteristics of DOM noise were determined 
using a minimum bias data set with pulses recorded by 120 



IceCube DOMs. The observed time difference between sequen- 
tial noise hits deviates from an exponential distribution expected 
for a Poissonian process (see Fig. The inset shows a hit 
sequence from a single DOM. Photomultiplier related after- 
pulses, wh ich occur wi t h » 6 % probability on time scales of 
0.3-1 1/is (Abbasi et al.L |2010|) . cannot explain the high occu- 
pancy bursts. We infer that the bursts are caused by an event 
within a DOM, but external to the electron amplification of the 
photomultiplier. 

Th ese observations are consistent with a study by iMeved 
who showed that a photomultiplier with bialkali pho- 
tocathode produces bursts that increase in rate and size as the 
photomultiplier is chilled. Such behavior could result from in- 
creased efficiency for radiative decay of excited states in the 
glass. In addition, Richardson's law describes the increase in 
thermal emission with photomultiplier temperature. The deploy- 
ment of IceCube DOMs on vertical strings places them in en - 
vironments ranging from -43 °C to -20 °C (IPrice et all 12002), 
warming with depth. DOM temperatures are some 10 °C warmer 
due to the energy dissipated in the electronics as monitored by 
a sensor mounted on the mainboard, but the photomultiplier and 
DOM glass temperature is somewhat uncertain. 

To confirm these temperature patterns, we divide the DOM 
noise into two contributions. The first is the rate of random ar- 
rivals as determined by fitting the slope of the interval distribu- 
tion as in Fig. [2] The second is the rate of events contributing 
to the excess of short intervals. These contributions are further 
divided into six temperature bands, and displayed in Fig. [3] The 
fitted excess cont ributions are then compared to an empirical ex- 
ponential ansatz (lMeveitl2010l) . 



r(T) = G-A c -. 



(3) 



where T is the absolute temperature, Ac = 0.055 m 2 is 

the cathode surface of the deployed R7081-02 photomultiplier 

dHamamatsuL 120071) and G = 5 x 10 4 /m 2 /s is a fixed constant 
i — £ Ia x . /fo7TTT3\ 



taken from iMeverf d2010l) . The fit results in T r * 115 K. The 
Poisson component is fitted to the Richardson-type law on ther- 
mal emission plus an ad-hoc constant noise term C 



r(T) = A-T 2 - e -*r+C 



(4) 



where k denotes Bolzmann's constant, W ~ 0.5 eV is the work 
function for Bialkali cathodes, and A, C are fit parameters. The 
curves match the observed temperature dependence of the noise 
rates fairly well, and support the hypothesis that the DOM noise 
is primarily due to a temperature dependent spectrum of bursts. 

The data acquisition was designed to reduce the noise rate by 
eliminating the excess hits, while keeping the random arrivals. 
The signal-to-noise ratio of the measurement can be improved by 
enforcing an artificial dead time r after every count, configured 
to 250 //s by a field programmable gate array in the DOM. This 
reduces the noise rate from 540 Hz to 286 Hz at the cost of some 
13% dead time for signal. The choice of 25 us optimizes sen- 
sitivity to the Lawrence-Livermore model (To tani et all 1 199% 
for distances up to 75 kpc, when neglecting the effect of after- 
pulses following the signal. A dead time r > 110/zs guarantees 
that the 4-bit counters do not overflow. The rate decrease due to 
dead time can be corrected for, however, the corresponding un- 
certainty increases once the measured rate approaches 1/t. An 
improved data acquistion that would avoid distortions of the rate 
measurement for supernova distances < 0(1 kpc) is under dis- 
cussion. Various dead time implementations, e.g. schemes mask- 
ing hits that arrive within the dead time caused by a previous hit 
or schemes allowing each hit to restart the dead time, were tested 
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ton induced compton electrons arising in the de-excitation of gi- 
ant O* resonances in the neutral current interactions vy + 16 O 
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Fig. 3. Measured average noise rates of 120 DOMs as function 
of DOM temperature. The excess of short time i nterval s due t o 
bursts (solid) is fitted to the empirical model of iMeveJ (12.010b . 
The Poissonian contribution (dashed) has not been corrected for 
the depth (and thus temperature) dependent contribution of at- 
mospheric muons (see Fig. [7). The dotted line is a comparison 
to the predicted thermal noise from Richardson's law plus a con- 
stant rate of C - 225 ± 6 Hz. 



v x + 16 O* -> v x + 15 0/ 15 N + n/p + y (Langa nke et all [l 996) 
have been included, we have not yet considered delayed f3 and 
j3y decays from excited nuclei. The cross section uncertainties of 
the reactions on protons and electrons are estimated in the refer- 
ences of Table[T]to be smaller than 1 %; uncertainties on oxygen 
reactions are hard to assess and the cross section may be only 
known up to a factor of two. 

Reactions producing electrons or positrons in the final state 
radiate N y Cherenkov photons along their flight path x, as long 
as their kinetic energies exceed the Cherenkov threshold of 
0.272 MeV. Integrating the Frank-Tamm formula between (300 
- 600) nm and accounting for the dispersion in ice, one obtains 
N y = (316 + 9) ctrT 1 • x (x in cm). For the inverse beta decay, the 
total average positron energy E e + is calculated from the average 
y e energy by the relation E e * ~ (£V, - 8) - - EyJ (Ey e + m p c 2 )] 

with 6 = (m^-m^-ml) /2m p , Due to the approximately quadratic 
energy dependence of the interaction cross section, the observed 
positron energies are on average higher than those of the incom- 
ing neutrino; the number of Cherenkov photons approximately 
rises with E\ (see Fig. |4j). 



with only slight differences observed. By eliminating the initial 
hits of the bursts, the noise rate can be reduced by up to 100 
Hz. Other optimizations may be possible but require a thorough 
understanding of the effect on signal hits. Design and implemen- 
tation of a new data acquisition system with more efficient noise 
rejection will be reported in a subsequent publication. 



4. Neutrino Interaction, Detection and Analysis 

In this section we first discuss the processes that lead to a de- 
tectable signal in IceCube, starting with the relevant neutrino in- 
teraction cross sections and continuing with the Cherenkov light 
production, propagation and detection. The real-time analysis in- 
troduced to monitor collective rate changes in IceCube's light 
sensors is described in the second part of this section. 



4. 1 . Cherenkov Photon Signal in IceCube 

Neutrinos of different species will be detected in IceCube via 
the interactions listed in Table Q] The table also includes the 
number of observed photon hits and the corresponding fractions 
as expected from the Garching model. The inverse beta reac- 
tion dominates in the ice. The small contribution by neutrino- 
electron scattering processes poses a challenge to the detection 
of the deleptonization peak. Note that the v e and v e cross sec- 
tions on lfi O are strongly energy dependent due to the high reac- 
tion thresholds of 15.4 MeV and 11.4 MeV, respectively. While 
their contribution to the hit rate in case of the Garching model 
is small (« 3 %), th e contribution can be a s high as 20 % for 
a 40 M progenitor dSumivoshi et all 120071) with average neu- 
trino energies of 25 MeV and beyond. The v e cross sections on 
the rare isotopes I8 0, 17 and on jH, with ls O giving the dom- 
inant contribut ion, add a small signal due to their low reaction 
thresholds (lHaxtonlll987l) . As the cross se ctions are only given 
for electron energies between (5-13) MeV (Haxton, 1999), they 
were extrapolated assuming an E\ dependence. While the energy 
deposition due to positron annihilation, neutron capture and pho- 
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Fig. 4. Example for an anti-electron neutrino energy spectrum 
with a = 2.92 pred icted for an 8.8 solar m ass O-Ne-Mg core 
supernovae collapse (Htidep ohl et all l2010) one second after the 
onset of the burst (solid line). Also shown is the cross section 
weighted energy spectrum of produced positrons and electrons 
from inverse beta decay (dashed line) as well as a measure of 
the detectable energy, which is proportional to the number of 
Cherenkov photons N y in the (300 - 600) nm range (dotted line). 
The relation N y - a - E e with a = 178 MeV -1 is used. 



The mean travel path for (9(10 MeV) electrons from v e and 
positrons from v e , including secondary leptons with energies 
above the Cherenkov threshold as well as positron annihilation, 
was determined with a GEANT-4 Monte Carlo simulation. The 
linear relationship x = (0.560 ± 0.005 (stat.) ± 0.030 (syst.)) cm 
■E e + 1 MeV was found for positrons. The corresponding relation- 
ship for electrons was determined to be consistent within errors. 

The optical scattering and absorption in glacial ice at the 
South Pole has been studied extensively (Acker mann et all 
120061) by the AMANDA and IceCube Collaboration with pulsed 
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Table 1. Major neutrino reactions. 



Reaction # Targets # Signal Hits Signal Fraction Reference 



v e + p — > e + + n 


6 


10 37 


134 k (157 k) 


93.8 % (94.4%) 


Strumia & Vissani (2003) 


v c + e~ — > v e + e~ 


3- 


10 38 


2.35 k (2.25 k) 


1.7% (1.4%) 


Marciano & Parsa (2003) 


v e + e~ — > v e + e~ 


3- 


10 38 


660 (720) 


0.5 % (0.4 %) 


Marciano & Parsa (2003) 


v^t + e~ -» + e" 


3- 


10 38 


700 (720) 


0.5 % (0.4 %) 


Marciano & Parsa (2003) 


Vfi+T + e~ -» v v+T + e" 


3' 


10 38 


600 (570) 


0.4 % (0.4 %) 


Marciano & Parsa (2003) 


v c + lfi O -> e" + X 


3' 


10 37 


2.15 k (1.50 k) 


1.5% (0.9%) 


Kolbe et al. (2002) 


v e + 16 -> e + + X 


3' 


10 37 


1.90 k (2.80 k) 


1.3% (1.7%) 


Kolbe et al. (2002) 


Van + 16 -> v M + X 


3- 


10 37 


430 (410) 


0.3 % (0.3 %) 


Kolbe et al. (2002) 


v e + ' 7/18 0/jH->e-+X 


6 


10 34 


270 (245) 


0.2% (0.2%) 


Haxton(1999) 



Notes. The approximate number of targets in a 1 km 3 ice detector, the detected number of hits at 10 kpc distance and their fraction in stars are given 
in the second, third and fourth column, respectively. In order to indicate the effect of neutrino oscillations in the star, signal hits and fractions are 
presented both assuming a normal neutrino hierarchy (Scenario A) and - in brackets - assu ming an inverted hierarchy (Scenario B). The numbers 
are taken from the Garching model using the equation of state by Latti mer & Swestvl (1991), integrating over 0.8 s and averaging over the neutrino 
incidence angle. 



and continuous light sources embedded in the ice with the neu- 
trino telescope. The detectors span depths ranging from (1300 
- 2500) m in the ice where the scattering coefficient varies by a 
factor of seven and absorptivity can va ry by a factor of three de- 
pending on the wavelength. The data ( Bra mall et al.L 120051) are 
consistent with the variations in dust impurity concentration seen 
in ice cores sampled at other Antarctic sites to track climatolog- 
ical changes. In the simulation applied for this paper, the ice is 
assumed to be homogeneous in the horizontal plane despite an 
observed slight tilt. 

We use two alternative procedures to calculate the num- 
ber of detected signal hits from the number of neutrinos 
crossing the detector: the first approach is based on separate 
simulations of particle interactions, Cherenkov photon creation, 
propa gation and detection, the second GEANT-3.21 GCALOR- 
based (Ze itnitz et all 1 1994b simulation combines all the steps in 
one program. 

IceCube's standard simulation of ph oton propagation withi n 
the ice relies on predetermined tables (Lundberg et al., 2007), 
created to track photons across the Antarctic ice. The tables store 
the detection probability and the arrival time distribution for 
given source and detector locations as well as their orientation. 
It includes the source wavelength, angular and intensity infor- 
mation, DOM parameters such as the glass and gel transparency 
and the quantum efficiency of the photomultiplier tubes. It also 
contains information about the ice such as the depth-dependent 
absorption and scattering lengths. 

The signal hit rate per DOM for a specific reaction and target 
is given by: 

DU\ "target ^SN^) f°° , v f°° , v 
K\t) - fdeadtime "= 

4nd 2 E y (t) Jo Jo 
X ^-(E e ,E v )N Y (E e )Vf f(E v ,E~ v ,a v ,t) , (5) 

where n tmget is the density of targets in ice, d is the distance of the 
supernova, L^ N (t) its luminosity, f(E v ,E v ,a v ,t) is the normal- 
ized neutrino energy distribution defined in Eq.|2]and E e denotes 
the energy of electrons or positrons emerging from the neutrino 
reaction. V£ ff denotes the effective volume for a single photon 
and N y (E) as 178 • E e is the energy dependent number of radi- 
ated Cherenkov photons; their numerical values depend on the 
selected wavelength range, chosen as (300 - 600) nm through- 
out this paper. The artificial dead time t (see Sect. \3.2\ reduces 



the total rate of hits. Comparing the observed signal, defined as 
the net increase over the nominal noise level, to the full rate of 
signal hits defines the dead time efficiency edeadtime- The approxi- 
mate expression edeadtime ~ 0.87/(1 + rsN • t) is found as function 
of signal rate rsN by adding Poissonian signal to the measured 
sequence of noise hits and applying a non-paralyzable dead time 
T = 250jt/s. 

The single photon effective volume varies strongly with the 
photon absorption. As a first approximation, V^ s can be esti- 
mated by the product of the Cherenkov spectrum and DOM sen- 
sitivity weighted absorption length (« 100 m), DOM geometric 
cross section (0.0856 m 2 ), Cherenkov spectrum weighted optical 
module sensitivity (« 0.071), average angular sensitivity includ- 
ing cable shadowing effects (^0.32), and the fraction of single 
photon hits passing the electronic DOM threshold (w 0.85). 

V^ ff was simulated by randomly placing 10 7 photons with 
(300 - 600) nm wavelengths within a sphere of radius 250 m 
around each DOM. We made the simplifying assumption that 
the Cherenkov light arrives at the DOMs isotropically from all 
angles. Note that the directions of positrons from the dominant 
inverse beta decay reaction are very weakly correlated with those 
of the incoming neutrinos. 

V?F was determined as function of depth in the ice (see 
Fig. |5). Averaging over all DOMs in one string one obtains 
Vf = 0.163 ± 0.004 (stat.) ± 0.020 (syst.)m 3 . The systematic 
uncertainty is discussed in Sect. 15.51 

The energy dependent effective volume Vf 1 for detecting 
an electron or positron is obtained by multiplying with the 
number N y (E) of Cherenkov photons. The mean number of pho- 
tons recorded by an optical module averaged over energy is then 

given by A^ etect = e deadtime • < teract • Vf> where < nteract is the 
neutrino density. For positrons with a cross section weighted av- 
erage energy of e.g. £e+=20 MeV (see Fig.|4|i one would obtain 

the average effective volume Vf = (29.0±3.8) m 3 /MeV • ~ 
(580 ± 80) m 3 for standard efficiency DOMs. This volume corre- 
sponds to an envisioned sphere of « 5.2 m radius centered at the 
optical module position, with full sensitivity inside and zero out- 
side. With 5 160 optical modules deployed, IceCube thus roughly 
corresponds to a dedicated 3.5 Mton supernova search detector 
in terms of geometry. Due to the presence of noise, a fair com- 
parison in terms of statistical accuracy needs to take into account 
the signal over background ratio as function of time and distance. 
To give an example, a study of the initial 380 ms of the burst in 
the Lawrence Livermore model (see Table [4]i at distances of 10 
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Fig. 5. Effective volume Vy S per DOM (left axis) for detection of 
Cherenkov photons with (300 - 600) nm wavelength plotted as a 
function of depth. The effective positron volume can be read off 
the right axis. DeepCore strings are not included in these plots. 



kpc (5 kpc) would require a 0.45 (1.6) Mton background free 
detector to statistically compete with IceCube. 

The second approach was to apply a GEANT-3.21 GCALOR 
based simulation of individual events that includes v e and v e 
on protons, electrons and 16 / 17 / 18 0, positron annihilation and 
neutron capture, the photon propagation in the ice including 
the effect of dust layer s, detector geometry, and the DOM 
response dRichardl 120081). The v e + P ~ > e + + n cross section 
parametrization of IVoeel & Beacom) (1 19991 ). which is in good 
agreement with IStrumia & Vissanil d2003l) . was used. Positron 
annihilation and hydrogen capture of neutrons produce photons 
of 0.5 1 Me V and 2.22 MeV energy, respectively. These add, pre- 
dominantly by Compton scattering and subsequent Cherenkov 
emission, w 1 MeV to the recorded energy. Rates from neutrino 
interactions on electrons reveal a 20% dependence on the 
incoming neutrino direction due to the small angle between 
neutrinos and scattered electrons and a directional dependence 
of the DOM efficiencies. Fig.|6]shows the clustering of detected 
inverse beta neutrino interactions at the position of the detector 
strings to visualize the effective volumes. The use of events with 
two or more DOMs detecting photons from the same positron to 
improve upon IceCube's sensitivity at large supernova distances 
and to track relative changes in the average neutrino energy will 
be discussed in a future paper. If several photons arrive close 
in time at the same DOM they will be counted as one hit; if 
one of the photons is delayed by scattering it will be rejected 
by the artificial dead time requirement. The two independent 
approaches for the determination of the detected number of 
events agree within 10 %. 

One may obtain a rate estimate from measured data by 
scaling the 11 events Kamiokande-II observed during the 
Supernova SN1987A neutrino burst to IceCube's effective 
v olume of Antarctic ice. A ssuming the v energy spectrum 
of lVissani & Paglioroli (2009), accounting for the Kamiokande- 
II energy threshold and positron detection efficiency, and tak- 
ing into account the loss due to IceCube's artificial dead time 
we determine a signal expectation of 113 + 36 detected photons 
per IceCube module within the first 15 s for a SN1987A like su- 
pernova near the galactic center at 10 k pc distance. The results 
are co nsistent with earlier simulations dFeseii 2004; Jacobsen, 
1996) performed for AMANDA that assumed homogeneous ice, 
after correcting for the different photomultiplier sensitive areas, 
optical module transparencies and dust layers in the ice. 
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Fig. 6. Detected neutrino inverse beta decay interaction vertices 
projected onto the horizontal plane based on a GEANT-3.21 sim- 
ulation with 10 million neutrino interactions. 



4.2. Real-Time Analysis Method 

The analysis monitors the collective rate increase Afi in all 
DOMs induced by Cherenkov photons uniformly distributed in 
the ice. As discussed in Sect. 14. II the photons are radiated by e* 
produced by reacting supernova neutrinos. Counting Nj pulses 
during a given time interval At, rates r-, = Nil At for DOM i, are 
derived. The index i ranges from 1 to the total number of oper- 
ational optical modules A^dom- With sufficiently large At's, the 
distributions of the r, 's can be described by lognormal distribu- 
tions that, for simplicity, are approximated by Gaussian distri- 
butions with rate expectation values (r,) and corresponding stan- 
dard deviation expectation values (<x,). These expectation values 
are computed from moving 300 s time intervals before and af- 
ter the investigated time interval. Shorter time intervals reduce 
the sensitivity of the analysis. At the beginning and the end of a 
SNDAQ-run, asymmetric intervals are used. The time windows 
exclude 30 s before and after the investigated bin in order to re- 
duce the impact of a wide signal on the mean rates. 

The most likely collective rate deviation Afi of all DOM 
noise rates r,- from their individual (r,)'s, assuming the null hy- 
pothesis of no signal, is obtained by maximizing the likelihood 



A'dom 



£(A/i) = n 



i 



1=1 



V2^V;> 



, (r,-«r,) + e,A / /))\ 

ex P ( ^r~a ) 



(6) 



Here e, denotes a correction for module and depth dependent 
detection probabilities. An analytic minimization of - In X. leads 
to 



a 2 v tiin-in)) 



with an approximate uncertainty of 



(No, 



^ i=i 



(7) 



(8) 



Note that A/i has the structure of a weighted average sum: each 
optical module contributes with the deviation of its expected 
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noise rate weighed by e,7(cr,) 2 . Assuming uncorrelated back- 
ground noise and a large number of contributing DOMs, the sig- 
nificance ^ = Ajt//cr AA , should approximately follow a Gaussian 
distribution with unit width centered at zero. In practice, the 
width turns out to be larger (see Sect. 15.31 ). The likelihood that a 
deviation is caused by an isotropic and homogeneous illumina- 
tion of the ice can be calculated from the ^-probability 

-2xo-4,-g( ^M±aM )'. < 9 , 

In order to suppress high rate deviations due to a temporary 
malfunction of individual detector modules, we reject supernova 
candidate events with a ^^-probability < 0.001. 

For short time bases Af, the Gaussian approximation is no 
longer valid and Poissonian probabilities must be used. The col- 
lective rate deviation can then be obtained from the equation 

= — ^ — = ^ 1 <^rr^7 - • (10) 

which can no longer be solved analytically. The same goes for 
the corresponding uncertainty, which is derived by identifying a 
drop in lnX by 0.5. The required numerical minimization pre- 
vents an online analysis of the raw data in 2 ms time intervals. 
However, fine time data in intervals of 30 s before and 60 s after 
a trigger are transmitted by satellite to perform a more detailed 
analysis offline. For instance, the onset of the neutrino emission 
can be determined with better than 5 ms accuracy for supernovae 
with less than 15 kpc distance conservatively assu ming low mass 
O-Ne-Mg core supernov ae ( Hude pohl et all l2010h. For simila r 
studies see lPapriaroli et al.l(l2009al) and lHalzen & Raffelil(l2009l) . 
This information can then be used to triangulate the supernova 
direction with other neutrino experiments. 

The optimal time base Af to detect faint signals depends on 
the expected signal shape. A simple generic description incor- 
porates a fast 0(1 ms) rise of the neutrino flux followed by an 
exponential decrease as expected during proto neutron star cool- 
ing with a time constant of r » 3 s. Maximizing 

Afi Jq e-rdt 

p= — (11) 

O" A/i VAf 

leads to Af » 1.26 r w 3.8 s. As the realtime analysis op- 
erates on bins of 0.5 s length, a time window length of 4 s 
has therefore been chosen as the best available setting for 
this pa rticular model a ssumption. Assuming the Livermore 
model (ITotani et al. , 1997), with a pronounced flux during the 
first seconds due to the high mass progenitor, the optimal time 
window is determined to be 1.6 s. 

To cover these model uncertainties, additional analyses with 
time bases of 0.5 s and 10 s are run in parallel with the one with 
Af = 4 s. The 0.5 s analysis aims at short neutrino bursts (i.e. 
from soft gamma ray repeater sources or from supernovae col- 
lapsing into a black hole). The 10 s time base accounts for the 
observed time window of the detection of the neutrinos from 
the supernova SN1987A by Kamiokande-II (Iffirataet al.lll987l) . 
By removing the cut on y 1 f° r the 0.5 s binning, the trigger 
has been made sensitive to partial illuminations of the detector. 
This gives the possibility to record hypothetical exotic particles 
emitting considerable amounts of light and thereby acting as a 
slowly moving source (such as ultra-heavy magnetic monopoles 
in some theories). 



The collective rate deviation A/u and its uncertainty cr^ in 
the time bases of 4 s and 10 s are calculated using sliding win- 
dows in 0.5 s steps and extracting the maximal significance. This 
procedure ensures that the signal detection efficiency is not re- 
duced by binning effects. 

5. Detector Performance 

In this section we will characterize the detector performance 
based on two years of data taking experience, discuss detector 
qualification criteria and summarize the expected systematic un- 
certainties. The data were taken with 22 operating strings (211 
days between Aug 2, 2007 and April 5, 2008) and with 40 op- 
erating strings (345 days between April 9, 2008 and April 15, 
2009). 

5.1. DOM Stability Requirements 

The stability of DOM noise rates is crucial for IceCube's sensi- 
tivity to detect supernovae. Faulty modules are removed from the 
analysis using automatic procedures that are applied in real time. 
In the 40 string configuration, 41 DOMs out of 2400 deployed 
DOMs showed no signal (« 1.7 %); all module rates fulfilled the 
requirement (r,) < 10000 Hz. Operational modules are removed 
from the analysis if they exhibit a variance (cr,) 2 much larger 
than the Poissonian expectation (r,) or high skewness In the 
very rare case, where the number of qualified modules drops be- 
low a threshold of 100, the corresponding time periods are dis- 
carded as a safeguard to prevent sending false alarms to SNEWS. 

The filter results in Table [2] show the excellent data quality 
of IceCube. Taking the 40 string configuration as an example, 
98 % of the disqualifications were due to just 1 1 DOMs. Thus 
the module disqualification has a negligible effect on the signal 
significance which changes as the square root of active DOMs. 

5.2. Long Term Stability 

The DOM rates are characterized by an exponential rate decrease 
over long time periods and a slight seasonal modulation. For the 
purpose of this analysis, the formula 

r(t) = r + c\e~* lT + c 2 sin(27r(f/year)) (12) 

represents the effects sufficiently well, as can be seen in Fig. [7] 
The decay of the rate is likely due to a decrease of tribolumines- 
cence in the ice with time, a byproduct of the freezing process. 
For DOMs that have been in the ice for more than 3 years, the 
fall time r exceeds 40 years, except for very deep DOMs where 
the freezing process takes longer (t « 4 years). In any case, the 
effect is negligible for the analysis which requires a stable rate 
within the analysis time window of lOmin. The slightly skewed 
rate distribution of a single DOM is better described by a log- 
normal distribution than by a Gaussian (see Fig. [8] with 250/is 
dead time applied). The average rates for standard efficiency and 
high efficiency DOMs are determined to be 286 Hz and 359 Hz, 
respectively. Thanks to the tight quality control, variations 
between DOMs of the same type are small showing standard 
deviations of 26 Hz and 36 Hz, respectively. The seasonal DOM 
rate modulation is assumed to arise fr om a change in the atmo- 
spheric muon flux ( Tila v et al.L l2009h . Fitting the time varying 
component, the parameter can be extracted (see Hig. [7J, which 
tracks the effect of dust layers similarly to what was observed 
in the determination of effective volumes (see Fig. |5). If one 
interprets the effect as being due to stratospheric temperature 
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Table 2. Module disqualification. 



Analysis 


Broadening cut 


Loss (%) 


Skewness cut 


Loss (%) 


Fraction of DOMs removed 


bin width (s) 










from analysis (%) 


0.5 


0.64 < (o-tfKr,) <4.0 


0.02 


|*| < 0.8 


0.01 


1.73 


4 


0.64 < (atfUn) <4.0 


0.05 


\s\ < 1.2 


0.02 


1.75 


10 


0.36 < (a-if/in) <6.25 


0.02 


14 < 1.7 


0.01 


1.73 



Notes. Disqualification requirements used in the online analysis and corresponding average percentage of DOMs that are rejected. The values 
were extracted from successful data taking runs with 40 IceCube strings covering the time from 2008/04/09 to 2009/04/19, which corresponds to 
an uptime of about 345 days. Most DOMs removed from the analysis are dysfunctional and provide zero rate. 
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Fig. 7. Top: Rate of a typical DOM as function of time cover- 
ing 556 days of lifetime as measured in 0.5 s bins (baseline sup- 
pressed). The line corresponds to a rate fit according to Eq.fTZl 
Bottom: Parameter q and estimated muon induced rate as func- 
tion of depth. The variation with depth is mostly due to the 
optical properties of the ice and muons ranging out. 



variations measured to modulate the muon flux by Ar w 8.3 % 
in 2008, the averaged muonic contribution to single DOM rates 
is c-il Ar m 12c2 ~ 16 Hz. As will be discussed in the following 
section, statistical fluctuations in the atmospheric muon rate, 
despite the small average muon contribution to the DOM rate of 
16 Hz, distort the significance spectrum considerably. 

5.3. Background Significance Distribution 

For purely uncorrected background and high statistics, the sig- 
nificance ^ is expected to be Gaussian distributed with width 
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Fig. 8. Rate distribution of a typical standard efficiency DOM 
taken over 29 consecutive days. Each measurement corresponds 
to 0.5 s integration time. Gaussian and lognormal fits are shown. 



cr = 1 (see Sect. I4.2l i. As can be seen in Fig. [9] the measured 
distribution is broader than expected and can be fairly well fitted 
by a Gaussian with width cr = 1.27. The broadening increases 
with the size of the detector and has reached cr - 1 .43 with 79 
operating strings. This broadening is due to non-Poissonian fluc- 
tuations in the number of hits deposited by atmospheric muons: 
highly energetic muons or muon bundles clustering in time leave 
correlated hits that will in general pass the x 2 cut. In the offline 
analysis, one can partly remove this effect as the number of muon 
induced coincident hits in neighboring DOMs is recorded for 
all triggered events. For the 79 (40) string configuruation, the 
broadening is thus reduced to cr = 1.06 (1.05). As this section 
describes the results of the online analysis, we do not apply this 
correction in the discussions below. 

An effective significance threshold of £ = 6.0 provides an 
internal trigger for testing the system one to two times per day, 
while a threshold at £ — 7.1 satisfies the SNEWS requirement 
of one false background trigger approximately once per 10 days. 
The Gaussian curve shown in Fig. [9] predicts one false back- 
ground trigger within ten years at a threshold at £ = 11. These 
thresholds are also depicted in Fig Q~2] The entries at £ = 8 and 
£, - 9.5 are due to test runs with artificial light sources. 

5.4. Future improvements 

Further optimizations may be applied to the data acquisition and 
analysis in the future, e.g. by incorporating a more sophisti- 
cated method to remove correlated noise, by excluding the bin- 
by-bin contribution of measured cosmic ray muon hits to the 
online rate measurement, by storing time stamps of all hits in 
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case of a significant alarm to e.g. improve o n the timing resolu- 
tion and to track the av erage neutrino energy (iBaum et all 1201 it 
iDemiroers et al.Ll201 lb . and by employing temporal templates in 
likelihood or cross-correlation studies. 




Significance ^ 

Fig. 9. Significance distribution in 0.5 s binning for a detector 
uptime of 556 days with 22 and 40 strings deployed. The two 
outliers at 8 and 9.5 occured during test runs employing arti- 
ficial light. The dashed line shows a Gaussian fit with <x = 1.27. 



5.5. Systematics 

There are three types of systematics relevant to this paper. The 
first type has to do with time dependent changes in the noise rate 
in all or a subset of DOMs that can mimic a supernova signal. 
These include high voltage variations, longterm trends such pho- 
tomultiplier aging, weather effects and other experimental effecs. 
For the supernova monitoring, realtime analysis as well as trig- 
gering, longterm trends are accounted for by calculating the rate 
expectation values (r,) and their standard deviation (cr,) from a 
rolling average of the noise rate for each DOM over 300 s on ei- 
ther side of the time of interest as was described in Sect. 14.21 ). 
The second type affects our understanding of the overall sensi- 
tivity of the detector. Ice properties, the wave length dependent 
quantum efficiency and the DOM thresholds all fall under this 
category. The third type is due to our current knowledge of rele- 
vant cross sections, the distance to supernovae, the neutrino-type 
dependent luminosity and energy, as well as oscillation effects in 
the star and in the Earth. 

Detector Stability and Environmental Effects 

The detector behaves very stably under normal operation. 
Periods during drilling, tests with artificial light sources and pe- 
riods with data acquisition problems as well as a few noisy mod- 
ules are excluded (see Sect. 15. U from the analysis. As discussed 
in 15.21 annual variations as well as shorter term modulations 
in the atmospheric muon flux change the observed rates which, 
however, are tracked by the rolling average. As hits from muons 
penetrating the detector are recorded simultaneously by the data 
acquisition system, they can be subtracted from the supernova 
rate measurements offline. Overall, the uncertainty on the super- 
nova sensitivity associated with the detector stability is estimated 
to be small (1.6%). 



The data were checked for other external sources of rate 
changes such as seismic activity and varying magnetic or electric 
fields as tracked by magnetometers and riometers at the South 
Pole. Only magnetic field variations show a slight, albeit in- 
significant, influence on the rate deviation of -1.3 ■ 10~ 6 Hz/nT. 
The influence is 30 times lower in IceCube than in AMANDA 
due to a wire mesh fi metal shielding in IceCube DOMs. 

Ice Properties and Sensitivity of the Detector 

Dust and air bubbles in the natural ice medium cause photons 
to scatter, while dust and the ice itself determine the absorp- 
tion lengt h. The range o f ice densities p lce = (919.6 ± 1.6) 
kg/m 3 dPrice et all 120021) reflects the 0.4 % density decrease 
due to the temperature increase between (1.4-2.4) km depth. 
Scattering dominates in the shallow ice above 1400 m and pos- 
sibly in the ice of the hole around the DOMs, which refreezes 
soon after deployment. Uncertainties in the optical properties of 
the hole ice are estimated to affect the effective volume deter- 
mination by < 1 %. More important are the uncertainties in the 
description of ice properties of the Antarctic glacier. 

The distributions of the photon arrival times and number 
of photons received at AMANDA modules from artificial light 
sources were used to derive the scattering length at different 
depths. Pulsed and continuous LED and laser sources give com- 
plementary measurements. The measurements of ice properties 
are consistent to within 6 % of each other including statistical 
and systematic uncertainties both for the scattering and absorp- 
tion measurements. The information on photon propagations is 
stored in tables contributing an estimated 1 % uncertainty due 
to the finite binning. Our knowledge of ice properties and cor- 
responding simulation methods continue to improve. Variation 
of DOM optical sensitivites and effects of the photomultiplier 
threshold on single photo-electron pulses lead to an « 10 % un- 
certainty. We assumed a (7 ± 3) % loss of light due to cables that 
shadow the photomultiplier surface. The effects discussed in this 
paragraph add up to an overall 12 % uncertainty. 

The track length of a positron or electron, including that of 
secondaries with kinetic energies above the Cherenkov threshold 
of 0.272 MeV, depends linearly on the initial lepton energy. The 
statistical uncertainty of the GEANT-4 calculation, including a 
systematic difference between electrons and positrons, is 0.3 %. 
The implementations of low energy electromagnetic processes 
have been cross checked between GEANT and NIST ESTAR- 
ICRU37 compilations. Go od agreement has be en found in par- 
ticular for electron ranges ( Ama ko et all |200"5|) . NIST quotes a 
(2 - 5) % systematic uncertainty on their implementation of elec- 
tromagnetic cross sections in the energy range relevant to super- 
novae. Event to event statistical fluctuations in the track length 
and in the number of Cherenkov photons (« 2 %) are negligible 
when investigating the ensemble of all DOMs. 

External Sources of Systematics 

The estimated uncertainties of the cross sections are listed in 
Table Q] Those associated with oxygen scattering processes are 
large and difficult to assess. Due to the strong energy depen- 
dence of I6 neutrino cross sections, the impact of this uncer- 
tainty depends on the energy spectra of particular models and 
the assumed oscillation scenarios. Processes involving only v e 
scattering are particularly affected. The total systematic uncer- 
tainty from detector effects and cross sections on the total rate of 
all neutrinos (electron neutrinos) is w 14 % (25 %). 
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Table 3. Summary of statistical and systematic uncertainties. 



Source of Systematics 


Uncertainties 




stat. (% ) 


sys. (%) 


Rate deviation in sliding average within 600 s buffer 


1.6 




Ice density as function of depth 




0.2 


Mean track length of electrons and positrons x of e* in ice 


-0.9 


5 


Uncertainty of the effective volume for photons 


1.3 


~ 12 


Uncertainty on the effect of artifical dead time on the signal 




3 


Uncertainty of the inverse beta-decay cross section 




< 1 


Uncertainty of the electron scattering cross sections 




< 1 


Uncertainty of the oxygen scattering cross section 




~ [3,20] 


Angle dependent MSW oscillation in the mantle and core of the Earth 




~ [0,-8] 



Notes. Estimated uncertainties on the total observed rate of all neutrino and anti-neutrino flavors. The earth effect strongly depends on the model 
and the neutrino angle of incidence. The uncertainty due to oxygen scattering depends on the neutrino energies and the neutrino channel. Assuming 
oxygen cross sections are known within a factor of two, rates may change between 3 % (Garching model) and 20 % (black hole model). For electron 
neutrino scattering alone, oxygen scattering accounts for about 40 % of the rate in the Garching model, with correspondingly large uncertainties. 



The distanc e to stars in our Galaxy is typically known to 
25 % accuracy dScheffler & Elsasseit Q998). However, the dis- 
tance of a supernova can in principle be measur ed by interpret- 
ing it s light curve with an accuracy of (5 - 10) % (Eas tman et all 
1996). Unfortunately, a considerable fraction of supernovae oc- 
curring in our Galaxy may be obscured by dust at optical wave- 
lengths. 

The uncertainties in the supernova collapse models are large 
and difficult to assess. The v e rate from the neutronization burst 
is largely independent of the progenitor mass; the corresponding 
uncertainties are estimated to be around 10 %; uncertainties aris- 
ing from neutrino o scillations are estimated to be below 5 % for 
a normal hierarchy (Kac helriess et all [2005). 

Oscillations in the Earth strongly depend on the incoming 
neutrino direction and may lead - depending on the neutrino 
hierarchy - to a maximal rate decrease of 8 % and 3 % during 
the cooling phases of the Lawrence Livermore and Garching 
models, respectively. The differences between various oscilla- 
tion scenarios may be as large as 30 % or even 50 % in the case 
of black hole formation. 



6. Performance Simulations 

In this section we will discuss the capability of IceCube to char- 
acterize details of the core collapse of massive stars and of the 
supernova remnant, as well as the insights IceCube may pro- 
vide into the properties of neutrinos and their interactions. There 
remain significant uncertainties in our understanding of the neu- 
trino emission from supernova explosions, necessitating com- 
parisons between several models to map the parameter space. In 
order to illustrate IceCube's performance, we will refer to spe- 
cific models chosen to span the possible range of supernova pro- 
genitor masses and neutrino energy spectra. We will also refer 
to more speculative models in order to demonstrate IceCube's 
high statistical precision in the detection of modulations of the 
neutrino light curve from astrophysical effects. 

When discussing the complete accretion and cooling phase 
exten ding to 15 s, we re fer to recent O-Ne-Mg core mod- 
els dHudepohl et all l2010h and the older Lawrence-Livermore 
model dTotani et all 1 19971) as examples with low and high 
progenitor star masses. The calculations consider only the 
radial dimension as a parameter. When discussing the first 
800 m illiseconds of the b urst we also refer to the Garching 
model dKitaura et all 120061) as an example for calculations with 



sophisticated transport mechanisms. Because it assumes only 
half of the initial star mass, (8 - 10) M instead of 20 M , it 
predicts fewer neutrinos than the Lawrence-Livermore model. 

In order to be compatible with other studies, we will usu- 
ally show experimentally predicted neutrino light curves for 
distances of 10 kpc, roughly corresponding to the center of 
our Galaxy. Depending on the m odel for the s upern ova precur- 
sor distribution, b etwee n 44% dAhlers et all l2009h and 53 % 
dBahcall & Piranl 1 19831) of all core collapse candidate stars 
in the Milky Way are expected to occur within this distance. 
Abou t 90 % of all s upern ovae are predi cted to occur within 15.4 
kpc dMirizzi et all 120061) to 17.5 kpc dAhlers et all l2009h dis- 
tance from the Earth. 

In the study of star matter oscillation effects, we restrict our- 
selves to the comparisons of the three scenarios A-C for neutrino 
hierarchy and 613 mixing angles that were introduced in Sect. [2] 
For some comparisons, we also show distributions with star mat- 
ter oscillations turned off. 

All simulations are performed for the final IceCube array 
with 4800 standard and 360 high efficiency DOMs. We assume 
that 2 % of the DOMs are excluded from the analysis, either be- 
cause they are not working or they give unstable rates. The back- 
ground noise was accounted for in two different way. For the de- 
termination of the significance and galaxy coverage, the simu- 
lated signal was randomized assuming a Poissonian distribution 
and added to noise data taken from experimental measurements 
and analyzed with the real-time reconstruction programs. For 
the simulation and comparision of various models, we added the 
calculated and randomized signal rates to the noise of the floor 
drawn from a Gaussian with mean value and standard deviation 
derived from data. 

Due to correlated pulses from radioactive decays and atmo- 
spheric muons, the measured sample standard deviation in data 
taken with 79 strings is a: 1.3 and * 1.7 times larger than the 
Poissonian expectation for 2 ms and 500 ms bins, respectively. 
It is possible to subtract roughly half of the hits introduced by 
atmospheric muons from the total noise rate in the offline anal- 
ysis, as the number of coincident hits in neighboring DOMs is 
recorded for all triggered events. We apply this correction to all 
Monte Carlo analyses described in this section, as this procedure 
lowers the standard deviation to (1.24- 1.32) sjj^i r,-, slightly de- 
pendent on the binning. 

Unless noted otherwise, we will use a likelihood ratio 
method to determine the range within which models can be 
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distinguished. From sets of several thousand test experiments, 
we will typically determine limits at the 90 % confidence level, 
while requiring that the tested scenario is detected in at least 
50 % of the cases. Note that the ranges obtained should be inter- 
preted as optimal as we assume that the model shapes are per- 
fectly known and only the overall flux is left to vary; we also 
disregard the possibility that multiple effects, such as matter in- 
duced neutrino oscillations and neutrino self-interactions, could 
co-exist and thus may be hard to disentangle. 

6.1. Expected Supernova Signal 

Evaluating Eq. [5] one obtains the rate spectra of Fig. [10] for a 
supernova at lOkpc distance. With a maximal signal-over-noise 
ratio of * 55 for the Lawrence-Livermore model, the neutrino 
burst can clearly be detected with IceCube. Also, the still hy- 
pothetical accretion phase lasting from (0 - 0.5) s can be sep- 
arated from the subsequent cooling phase with high statistical 
precision. The study of the cooling phase is limited by the pho- 
tomultipli er noise in particular fo r the case of the light O-Ne-Mg 
model by Hiidepohl et alj d2010t) . 
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Fig. 10. Expected rate distribution at lOkpc distance for the 
La wrence-Livermore mo del (dashed line) and O-Ne-Mg model 
by Hudep ohl et al.l d2010i) with the full set of neutrino opacities 
(solid line). . The 1 cr-band corresponding to measured detector 
noise (hatched area) has a width of about ± 330 counts. 



The oscillation scenario B for an inverted neutrino mass hi- 
erarchy shows the largest signal for the Lawrence-Livermore 
and Garching models because energetic v x will oscillate into v e , 
harden their spectrum and thus increase the detection probabil- 
ity. The scenario without any oscillation is presented as a ref- 
erence and leads to the weakest signal. Scenario A (normal hi- 
erarchy) and Scenario C (very small #13 < 0.09°) are hard to 
distinguish due to their very similar effect on neutrino mixing. 

Clear differences between the oscillation scenarios in abso- 
lute rate and shape appear in Fig. [TT] Assuming that the model 
shapes are known but not necessarily the overall normalization, 
the inverted hierarchy can be distinguished from the null hypoth- 
esis of a normal hierarchy up to distances of 16 kpc. 

6.2. Significance and Galaxy Coverage 

The simulation of an expected signal from a supernova within 
the Milky Way has to take into account the number of likely pro- 



genitor stars in the Galaxy as a function of the distance from 
Earth. The expected significances of supernova signals accord- 
ing to the Lawrence-Livermore model for three oscillation sce- 
narios are shown in Fig. Q~2] For this particular model, the sig- 
nificances for the 4 s and 10 s binning turn out to be approx- 
imately 20 % and 50 % lower than for 0.5 s, respectively. For 
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Fig. 12. Significance versus distance assuming the Lawrence- 
Livermore model. The significances are increased by neutrino 
oscillations in the star by typically 15 % in case of a normal hi- 
erarchy (Scenario A) and 40 % in case of an inverted hierarchy 
(Scenario B). The Magellanic Clouds as well as center and edge 
of the Milky Way are marked. The density of the data points 
reflect the star distribution. 



the graph, the superno va progenitor distribution predicted by 
Bahcall & Piranl l l 19831) was used. For the Magellanic Clouds, 
which contain roughly 5 % of the stars in the Milky Way, a uni- 
form star distribution along the diameters of the galaxies was 
assumed for simplicity. 

IceCube is able to detect supernovae residing in the Large 
Magellanic Cloud (LMC) with an average significance of (5.7 ± 
1 .5) <x in a 0.5 s binning, assuming the Lawrence-Livermore 
model. The uncertainty reflects different oscillation scenarios. 
Supernovae in the Small Magellanic Cloud (SMC) can be de- 
tected with an average significance of (3.2 + 1.1) <x and will in 
general not trigger sending an alarm to SNEWS, as indicated 
by a horizontal line in Fig. [12] IceCube will observe supernovae 
in the entire Milky Way with at least a significance of 12 <x at 
30 kpc distance. 

6.3. Onset of Neutrino Production 

The analysis of the deleptonization peak that immediately fol- 
lows the collapse is of considerable interest, since its magnitude 
and time profile are rather independent of the initial star mass 
and of the nuclear e quation of state; the variation is estimated 
by dKeil etalll2003h to be around 6 %. Thus the electron neu- 
trino luminosity may be used as a standard candle to measure 
the distance to the supernova. 

As the deleptonization peak lasts for only 10 ms, the data are 
evaluated in the finest available time binning of 2 ms, as depicted 
in Fig. [TT] The deleptonization signal is detected by the elastic 
v e + e~ — > v e + e~ reaction with a cross section times the number 
of targets ss 50 times smaller than for the v e + p — > e + + n in- 
teraction. As the v e flux rises rapidly following the collapse, the 
deleptonization peak remains almost completely hidden, espe- 
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Fig. 11. Top: Expected rate distribution at lOkpc supernova distance for oscillation scenarios A (normal hierarchy) and B (inverted 
hierarchy). Fluxes and energies in th e left plot are taken f r om th e Lawrence-Livermore model and in the right plot from the Garching 
model using the equation of state of Lattim er"& Swestvl d 199 lh . Scenario C (not shown) is almost indistinguishable from Scenario 
A. The case of no oscillation is given as a reference. Bottom: Expected average signal rate distribution at lOkpc distance in finest 
2 ms binning for Scenarios A and B using the Garching model; the unlikely case of no oscillation is given as a reference. The 
left plot shows the expected v e induced signal. As can be seen from the right plot, the signal is no longer apparent, once the large 
contribution due to the inverse beta decay and the expected DOM noise are added. The 1 cr-bands corresponding to measured 
detector noise (hatched area) have a width of about ±215 counts for a 20 ms binning and + 70 counts for a 2 ms binning. 



cially when neutrinos oscillate in the star. In this case the subtle 
structure may be resolved only for distances d < 2 kpc. 

Largely independently of the model, each oscillation sce- 
nario shows a characteristic slope of the rate increase around 
the deleptonization peak. Quantifying this by a series of several 
thousand simulations for the Garching and Lawrence-Livermore 
models and considering oscillation Scenarios A-C and the case 
of no oscillation, it is possible to establish the inverse hierar- 
chy (Scenario B) w.r.t. the normal hierarchy (Scenario A) with 
90% C.L. for distances d < 6 kpc (corresponding to 21 % of 
all pr ogenitor stars, wh en accounting for the effect of spiral 
arms dAhlers et alll2009h 1. 

6.4. Shock Waves 

For an inverted hierarchy (Scenario B), the rate distribution 
should reveal the effects of forward and backward moving shock 
waves traveling through the collapsing star during the cooling 
p hase, (3 - 10) s aft er bounce. Assuming the specific model 
of Tomas et alJ d2004l) (see Fig. [TBI, scenarios with a static den- 



sity profile and one forward shock wave can be distinguished at 
90% C.L. up to distances of 13 kpc; the distance reduces to 10 
kpc in a scenario with one forward and one reverse shock wave 
(not shown). 

6.5. Quark Star and Black Hole Formation 

IceCube is particularly well suited to study fine details of the 
neutrino flux as function of time. As an example, Fig. [Pfl shows 
a simulation based on the prediction of iDasgupta et al.l {2010) 
for the formation of a quark star. The model predicts a sudden 
spike in the v e flux lasting for a few ms while the neutron star 
turns to a quark star; the time of the QCD phase transition can be 
determined with sub-ms accuracy. The likelihood ratio test gives 
a deviation larger than 5 <x from the hypothesis of no quark star 
formation for distances up to 30 kpc. Height and shape of the 
peak depend on the neutrino hierarchy. Scenarios A and B can 
be distinguished at 90 % C.L. up to distances of 30 kpc. 
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Fig. 13. Effect of a forward moving shock wave applied to super- 
nova at 10 kpc distance, modelled according to the Lawrence- 
Livermore model assuming an inverted hierarchy and 613 > 0.9° 
was assumed. A forward shock wave can be distinguished from 
a static density profile and the case of no star matter effect. The 
1 cr-band corresponding to measured detector noise (hatched 
area) has a width of about ± 1 150 counts. 



o 



7000 



6000- 



5000- 



4000 



3000 



No Oscillations 
■ Scenario A (NH) 
Scenario B (IH) 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
Time Post-Bounce [s] 



0.45 



Fig. 14. Comparison of the neutrino light curve with quark- 
hadron phase transition for a 10 M Q progenitor at 10 kpc dis- 
tance. Three neutrino oscillation scenarios are shown (see leg- 
end). The observation of the sharp v e induced burst, 257 ms 
< t < 261 ms after the onset of neutrino emission, would con- 
stitute direct evidence of quark matter. The hatched 1 cr-band 
corresponding to detector noise has a width of about ± 70 counts. 



Fig. [T31 show s a s imulation based on the prediction 
of Sumivoshi et al. (2007) for the formation of a black hole fol- 
lowing a collapse of a 40 solar mass progenitor star. Neutrinos 
reach energies up to 27 MeV (v e and v e ) and 40 MeV (v^ and 
v T ), carry a correspondingly large detection probability and thus 
produce very clear evidence for the formation of the black hole 
after 1.3 s, when the neutrino emission is expected to fade ex- 
ponentially (not realized in th e simulation). For Fig. Q3J a hard 
equation of state (IShenL fl998) was chosen, leading to black hole 
formation after 1.3 s. This corresponding drop can be identi- 
fied at higher than 90 % C.L. for all stars in our Galaxy and the 
Magellanic Clouds. 
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Fig. 15. Expected neutrino signal from the gravitational collapse 
of a non rotating massive star of 40 solar masses into a black 
hole at 10 kpc distance for a hard equation of state dShenl [l998) 
following Sumivos hi et al.l ((2007). The 1 cr-band corresponding 
to detector noise (hatched area) has a width of about + 70 counts. 



6.6. Neutrino Hierarchy Sensitivity and Rate Summary 

The number of standard deviation with which normal and 
inverted v hierarchies (Scenarios A and B) can be distinguished 
are plotted in Fig. [16] as function of the supernova distance 
for selected models. The values represent the optimal cases 
when model shapes (but not necessarily the absolute fluxes) are 
perfectly known. Table [4] lists the number of neutrino induced 
photon hits that would be recorded by IceCube on top of the 
DOM noise for various supernova models. Note that the number 
of expected signal hits scales with 1/distance 2 ; the dependence 
of the detection significance as function of distance can be read 
from Fig. [12] 



100 



10 



2? 
x 





L 


1 


1 1 






















































is 












% 














































5o 


s 






























h £>, L-i 










J.'"; 
■o — 












II 


--- 1 - 











10 15 
Distance [kpc] 



20 



25 



30 



Fig. 16. Number of standard deviation with which scenarios A 
(normal hierarchy) and B (inverted hierarchy) can be distin- 
guished in at least 50 % of all cases as function of supernova 
distance for some of the models listed in Table [4] A likelihood 
ratio method was used assuming known model shapes. 
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Table 4. Expected rates. 



Model 


Reference 


Progenitor 


#v's 


#v's 






mass (M Q ) 


t < 380 ms 


all times 


"Livermore' 1 


(Totani et al., 1997) 


20 


0.174 x 10 6 


0.79 x 10 6 


"Garching LS-EOS Id" 


(Kitaura et al.. 2006) 


8-10 


0.069 x 10 6 




"Garching WH-EOS Id" 


(Kitaura et al.. 2006) 


8-10 


0.078 x 10 6 




"Garching SASI 2d" 


(Marek et al.. 2009) 


15 


0.106 x 10 6 




"1987Aatl0kpc" 


(Pasliaroli et al.. 2009b) 


15-20 




(0.57 ± 0.18) x 10 6 


"O-Ne-Mg Id" 


(Hiidepohl et al.. 2010) 


8.8 


0.054 x 10 6 


0.17 x 10 6 


"Quark Star (full opacities)" 


(Daseupta et al.. 2010) 


10 


0.067 x 10 6 




"Black Hole LS-EOS" 


(Sumivoshi et al.. 2007) 


40 


0.395 x 10 6 


1.03 x 10 6 


"Black Hole SH-EOS" 


(Sumivoshi et al., 2007) 


40 


0.335 x 10 6 


3.40 x 10 6 



Notes. Number of recorded DOM hits in IceCube (as #v's) for various models of the supernova collapse and progenitor masses assuming a distance 
of 10 kpc, approximately corresponding to the center of our Galaxy. A normal neutrino hierarchy is assumed. 



7. Conclusion 

A high statistics observation of the supernova neutrino flux 
would provide valuable information on astrophysics and the 
properties of neutrinos. IceCube was completed in December 
2010 and monitors m 1 km 3 of deep Antarctic ice for particle 
induced photons with 5160 photomultiplier tubes. Since 2009 it 
supersedes AMANDA in the SNEWS network. With a 250 /is ar- 
tificial dead time setting, the average DOM noise rate is 286 Hz. 
The rates remain constant over time with a small modulation in- 
duced by changes in the atmospheric muon flux; they hardly vary 
across the detector once the DOMs have been frozen in for a 
sufficiently long period. The data taking is very reliable and cov- 
ers the whole calendar year, including periods when new strings 
were deployed. The uptime has continuously improved toward 
a goal of > 98 % and reached 96.7 % in 2009. IceCube's sensi- 
tivity corresponds to a megaton scale detector for galactic super- 
novae, triggering on supernovae with about 200, 20, and 6 stan- 
dard deviations at the galactic center (10 kpc), the galactic edge 
(30 kpc), and the Large Magellanic Cloud (50 kpc). IceCube 
cannot determine the type, energy, and direction of individual 
neutrinos and the signal is extracted statistically from rates that 
include a noise pedestal. On the other hand, IceCube is currently 
the world's best detector for establishing subtle features in the 
temporal development of the neutrino flux. The statistical uncer- 
tainties at 10 kpc distance in 20 ms bins around the signal max- 
imum are about 1.5 % and 3 % for the Lawrence Livermore and 
Garching models, respectively. 

Depending on the model, in particular the progenitor star 
mass, the assumed neutrino hierarchy and neutrino mixing, the 
total number of recorded neutrino induced photons from a burst 
10 kpc away ranges between « 0.17 x 10 6 (8.8 M O-Ne-Mg 
core), * 0.8 x 10 6 (20 M G iron core) to * 3.4 x 10 6 for a 40 
M progenitor turning into a black hole. For a supernova in the 
center of our Galaxy, IceCube's large statistics would allow for 
a clear distinction between the accretion and cooling phases, an 
estimation of the progenitor mass from the shape of the neu- 
trino light curve, and for the observation of short term modula- 
tion due to turbulent phenomena or forward and reverse shocks 
during the cooling phase. The deleptonization peak associated 
with the neutron star formation, however, may be hard to ob- 
serve since the electron neutrino cross section in ice is small. 
IceCube will be able to distinguish inverted and normal hier- 
archies for the Garching, Lawrence-Livermore and black hole 
models for a large fraction of supernova bursts in our Galaxy 
provided that the model shapes are known and 9u > 0.9°. The 
slope of the rising neutrino flux following the collapse can be 



used to distinguish both hierarchies in a less model dependent 
way for distances up to 6 kpc at 90 % C.L. As in the case of the 
inverted hierarchy, coherent neutrino oscillation will enhance the 
detectable flux considerably. A strikingly sharp spike in the v e 
flux, detectable by IceCube for all stars within the Milky Way, 
would provide a clear proof of the transition for neutron to a 
quark star as would be the sudden drop of the neutrino flux in 
case of a black hole formation. 
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